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À computer program was written for use in sizing gas turbine counter- 
flow regenerators. The program simulates the semi-open gas turbine re- 
generative cycle under constant load conditions. Inputs to the program 
include compression ratio, turbine inlet temperature, component efficien- 
cies, horsepower output and regenerator core geometry. External re- 
generator core dimensions are varied and various outputs are computed 
such as thermal efficiency, mass flow rate, total relative pressure 
drop, regenerator effectiveness and temperatures and pressures at all 


key points in the cycle. 
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r. Introduction. 

The design of a heat exchanger involves a consideration of both 
the heat transfer rates between the fluids and the mechanical pumping 
power expended to overcome fluid friction and move the fluids through 
the heat exchanger. Kays and London state that for a heat exchanger 
operating with high density fluids, the friction power expenditure is 
generally small relative to the heat transfer rate, with the result 
that friction power expenditure is seldom of controlling influence (1). 
However, for low-density fluids, such as the gases flowing through a 
gas turbine regenerator, it is very easy to expend as much аесһап al 
energy in overcoming friction as is transferred as heat. The regenera- 
tor designer must therefore establish the proper balance between heat 
transfer rate and pressure drop. 

Two advantages of the gas turbine are compactness and low weight 
whereas heat exchangers are inherently rather bulky and heavy. Weight 
and size considerations must therefore be taken into account along with 
heat transfer rate and pressure drop balance if regenerators are to 
serve a useful purpose. Since the prime reason for regenerators is 
the improvement of overall thermal efficiency, it appeared that a fruit- 
ful approach to regenerator sizing and core selection would be considera- 
tion of the gas turbine cycle as a whole rather than the regenerator 
separately. Adjustments in internal regenerator core geometry and ex- 
ternal dimensions are then reflected directly in terms of changes in 
thermal efficiency rather than indirectly in terms of changes in pres- 
sure drop and heat transfer rate. A computer program was written for 
this purpose which essentially simulates the semi-open gas turbine re- 
generative cycle under constant load conditions. Inputs to the program 


are compression ratio, turbine inlet temperature, component efficiencies, 
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horsepower output, regenerator core geometry, etc. External regenerator 
core dimensions are then varied and various outputs are computed such as 
thermal efficiency, mass flow rate, total relative pressure drop, re- 
generator effectiveness and temperatures and pressures at all key points 
in the cycle. 

Since the program is based on fundamental relationships from thermo- 
dynamics, heat transfer and fluid mechanics, sections 2, 4 and 5 are a 


review of the basic theory utilized. 
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2. ‘Thermodynamics of Semi-Open Regenerative Cycle 


REGENERATOR 





OMPRESSOR == же 


TEMPERATURE, ENTHALPY 


ENTROPY 


FIGURE 1. SEMI-OPEN REGENERATIVE CYCLE 
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2.1 Description of the Actual and Ideal Cycles. 

Both the ideal and actual regenerative cycles are shown in Figure 
l. The state points indicating the ideal cycle are subscripted "5"; 
the nonsubscripted state points indicate the actual cycle. 

First describing the ideal cycle, air continuously enters the com- 
pressor in condition "Ils". It is compressed isentropically to state "2s" 
then enters the cold side of the regenerator, Theoretically, if the re- 
generator is large enough and the flow is slow enough, the air from the 
compressor can be heated reversibly to temperature "5s" at state "35", 
It next enters the combustor, some of it going around the outside of 
the combustor chamber proper and some furnishing oxygen for burning fuel 
which is continuously injected. Because of the temperature rise, the 
gases expand and enter the turbine in state '4s". After isentropic ex- 
pansion through the turbine, the hot exhaust gases enter the hot side of 
the regenerator. Again, if the regenerator is large enough and the flow 
slow enough the gases are cooled reversibly to temperature ''2s'' at state 
"6s" and exhausted. 


The ideal cycle can be described as follows: 


2:02.20.) H.-B. С 


5. =.) С. = on (2.2) 
Tu Ж ) ds = | (2.3) 


In the actual cycle, air continuously enters the compressor as in 


the ideal cycle; condition "1" is identical to condition "ls". 


ЕЕ й E 
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(255) 


It is then compressed in an irreversible adiabatic process to state '2" 
arriving at the same pressure as in the ideal case but with increased 


entropy due to fluid friction. 
R- Rs (2.6) 


There will of course be a pressure drop as the air flows through the 
regenerator and the air is heated only to the temperature at state "3" 
not "3s" because less heat is actually transferred than could be trans- 
ferred reversibly. The air and now gases undergo another pressure drop 
in flowing through the combustor; however they arrive at the turbine 


inlet at the same temperature as in the ideal case. 


les =e (2.7) 


To continue the analysis it is necessary to look ahead to state point 
"6", The pressure at this point can be considered constant because it 
is assumed that the turbine is exhausting into a large volume, most 


likely the atmosphere. 
R- Rs (2.8) 


There will of course be a pressure drop through the hot side of the re- 


5 must be higher than Pss if P. is to 


remain the same as in the ideal case. The gases will also experience 


generator which indicates that P 


another increase in entropy as they pass through the turbine. Finally, 
since less heat is actually transferred from the hot side of the re- 


generator than could be transferred reversibly, T6 is higher than Tos 


17 


2.2 Compressor. 
If the compressor ís considered a steady flow machine with kinetic 


energy change equal to zero, the works of the fluid for adiabatic flow 


аге: 


ideal: Wes = his- has (2.9) 
actual: We = Ена (2.10) 


The foregoing equations may be obtained from the steady flow energy 


equation (2). The compressor efficiency is then: 


J. = Ель иок _ his- Ars e 
C ACTUAL FLUID WORK. ^ ALL AS о (2.11) 


2.3 Combustion Process. 


А 


ЖТ (1 еотв) i/a gk 


Figure 2 Епегру Diagram for Combustor 
The energy diagram shown above for a steady flow combustor, is on 
the basis of one lb, air entering and aA lb, fuel/lb. air. The 
entering molecular internal energy and flow work of the air is h, and 
of the combustible is ИА А . The stored chemical energy in the fuel 
qgis measured from the zero datum of enthalpy. Departing energy is the 
enthalpy of (1+ Y¢/a )1b, of products and that part of the chemical 


energy (1- N come ) Vefa 3А that is not converted into molecular 


form. Combustor efficiency as used here is defined as: 


„ actual energy taken on by the gases 
Comb 


energy released during complete combustion 
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The resulting equation is 


hs + VGA (qos) (Aa e(l Mens) ua q 2.12) 
Then 


pe па- Аз (2.13) 
s/a (As + Heonsgr—ha 


The lower heating value of the fuel qg із used in gas turbine analysis 
as it is impractical to expand the gas to a point where water vapor from 
the products of combustion will condense and contribute any energy to 
the cycle (3). 

The enthalpy of liquid octane witn respect to a datum of octane 


vapor having h=0 at absolute zero is given by 


ns EMI ae BrU/Z.. FUEL (2.14) 


The enthalpy of other hydrocarbon combustibles can safely be taken as 


the same (the amount of energy is small) (1). 


2.4 Turbine. 
If the turbine is considered a steady flow machine with kinetic 
energy change equal to zero, the works of the fluid for adiabatic flow 


ideal: We = (+ (ще Асу) (2.15) 
actual: МА e (I4 sÁ) ha hs) (2.16) 


As for the compressor, the foregoing equations can be obtained from the 


steady flow energy equation. The turbine efficiency is then 


Те--Астия ELUID WORK _ has Аг CRUS 
EAL RK р ° 
10 Wo h 4s- Ах 
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2.5 Overall Thermal Efficiency. 
The overall thermal efficiency is defined as net work out divided 


by energy chargeable against the cycle. 
Nen = MO. 
се» (2.18) 


The energy chargeable against the cycle is the energy released by the 


complete combustion of the fuel. 
Ес Vr. $4 (2.19) 


The equation for thermal efficiency then follows from equations (2.10), 


(2.16), (2.18), and (2.19) 


Жек= ax Gela hc] we (2.20) 
£ /A. 44 


2.6 Horsepower Output. 
The equation for horsepower output follows directly from equations 


(2.18) and (2.19) 


БР- (а а/46 (4/4 (2.21) 


where m has units of lb. air per hour and 2544 BTU = one horsepower 


hour. 


2.7 Specific Fuel Consumption. 

Specific fuel consumption is defined as lb a fuel burned per unit 
work output. If the horsepower-hour is used as the unit of work then 
specific fuel consumption has the units of lb a fuel per horsepower-hour. 


It follows that 


SFC Ud (2.22) 
€ 4 
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S Regenerator Surface and Core Geometry. 

Certain geometrical relations are necessary in the application of 
the basic heat transfer and flow-friction data to a design problem. One 
particular set and form of these relations that has proven convenient 


for use in the computer program is given here. 


PANN e 
ЖАРА 
ADAAKAAA 1 

br ЛААЛ 

COLD ЖАЖА Ë 

HOT ДАДАЛ 


COLD 
HOT -^ 
COLD 


Figure 3 Simple Plate-Fin Counterflow Regenerator Core 
The following part of basic data are given for all surfaces in Kays 
and London (4). 
b = plate spacing 
Yn 2 flow passage hydraulic radius 


= ratio of total transfer area of one side of 
exchanger to volume between plates of that side 


a = plate thickness 
In addition, the ratio of fin area to total transfer area is normally 
given. 


The following geometrical factors are required for each of the two 
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sides, hot and cold, of the complete heat exchanger core. 


А = 
А zz 


Afr = 


y = 
S 


The equations 


total transfer area of one side of exchanger 
free flow area of one side 

frontal area of one side 

flow length of one side 

total exchanger core volume 


ratio of total transfer area of one side of exchanger 
to total exchanger volume 


ratio of free-flow to frontal area of one side of 
exchanger 


solid core volume 


below give the relations between surface and core 


factors for one side of the exchanger. Subscript "1" refers to any one 


side and "2" refers to the other side; factors without a subscript are 


common to both sides. 


Q 


r A (3.1) 


ү (3.3) 
ДД, ута 


ТҮЙІП m) 5 


2 4. l, ri 
hy rh +A З 


= (Аз) = A c Аһ) [Ас ) (3.6) 
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A, 4 (A L Aer), 


V © L pn duces = (Ac А.) 


- - — © Gee. 
mE 





4. Counterflow Exchanger Heat Transfer Theory 
4.1 Conduction Heat Transfer. 

When a temperature gradient exists in a body there is an energy 
transfer from the high temperature region to the low temperature region. 
The energy is transferred by conduction and the heat transfer rate per 


unit area is proportional to the normal temperature gradient, 
е Я 
x (&.1) 
When the proportionality constant із inserted, 


DEM 
4- —4AS- (4.2) 


where "k" is called the thermal conductivity of the material. 


i dcn 





хх |< 
Figure 4 Elemental Volume for One-Dimensional Heat 
Conduction Analysis 


If the steady state one-dimensional case is considered 


$7 -4A 4X (4.3) 
{= LA еы - (4.4) 


which integrates to give 





Figure 5 Sketch Illustrating Convection from a Plate 
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The convective heat transfer rate is related to the overall temperature 


difference between the wall and fluid and the surface area as follows: 


+= h AC (Tw - To) (4.5) 


where "h" is called the convection heat transfer coefficient. It may be 
anticipated that convection heat transfer will have a dependence on the 
thermophysical properties of the fluid (thermal conductivity, specific 
heat, density, viscosity), the speed at which the fluid is moving and 
some characteristic dimension. If a dimensional analysis is performed 
using h; u, ср, 2 ° ke» M and L as variables, the variables may be 
arranged into three independent dimensionless parameters. The para- 
meters depend on how the analysis was performed. One very useful set of 


parameters is 


^ с LUL =f, ee Da 7. (4.6) 
) 


WC, 


À function of the form 
P (Sz tv, R г) = 0 (4.7) 


must therefore exist. Such a function may be determined experimentally 


in the form 


Se Re 5 (Re) m 


for various heat transfer surfaces and geometries. The generalized 
COP NS 

grouping Sz iv is defined as the Colburn j-factor and the char- 

acteristic dimension "L" in "Re" is designated as the hydraulic dia- 


meter, "Dn. 
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A large number of curves expressing the relationship given by equation 


ome Fu Du _ C D. (6.9) 


(4.8) are published in Kays and London (1). One of the more common ex- 
perimental methods for obtaining these characteristics is described in 
a paper by Pucci, Howard and Piersall (4). 

After determining the Colburn j factor, St is easily found pro- 


viding the fluid properties necessary to determine Pr are obtainable. 


we ғ,” (4.10) 


The convection heat transfer coefficient is then evaluated by using 


au 
h= Jz Gy E 


4.3 Fin Effectiveness. 





Figure 6 Sketch Illustrating One-Dimensional Conduction 
Through a Rectangular Fin with Convection at 
the Surface 
Consider the one dimensional fin exposed to a surrounding fluid at 
a temperature Too . The temperature of the base of the fin is T ° 
Making an energy balance on an element of the fin of thickness dx as 


shown in figure: 


(energy in left face)=(energy out right face)+(energy lost by convection) 
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energy in left face - 4x A AX 


y 


energy out right face - ТЕТІ = АЗ 
x tx 


Let P equal the perimeter of the fin then: 
energy lost by convection = Age = А P kx (Т-Т,) 


Combining quantities, the energy balance yields 





- 4A (£1 + A 


AT T. ФЕ (т- ІК \=0 (4.12) 
Ах” 
Let GO =тТ- Т . Then (4.12) becomes 
и 
ek АР, 20 (4.13) 
Ax* АА 
One boundary condition is 
0- 8 0-1. 5X X*0 
For the case of a fin with an insulated end, dT/dx » 0 and dO /dx = 
Oat x = Z. 
Let 
i. (4.14) 
SA 
The general solution for (4.13) may be written 
онн АТ (4.15) 
Solving for constants C1 and с, the solution is obtained as 
= mA A^ X 
См: €———— — (4.16) 


The total heat lost by the fin must equal the energy conducted through 


the base of the fin 


$7-4^Íl 


— 
= 


27 


Then 


i =í A 6 ( | Lau pe а bes “> \ (4. 17) 
I+ ¢ к 
zl P4 A @ tak (me) 


Fin efficiency is defined as 
= actual heat transferred 
heat which would be transferred if entire fin area were at base 


temperature 


which becomes 


ТАРДА O, Z< (aa) boc) (4.18) 


pap 


To apply equation (4.18) in the computer program it is necessary to 


put it in terms of the heat exchanger core geometry. Referring to Figure 


5, 
Ә-21--25%%0 
A= sb 
so that 
24 
AS 


It is seen from Figure 3 that the fins of the heat exchanger are 
symmetrical with respect to a plane midway between the plates separat- 
ing the fluids. Therefore, the temperature gradient in the fin at this 
plane is zero and the above insulated end solution applies. 


Noting that Z ~ b/2, the fin efficiency becomes 


= _ 4 (124/45 AL (nas %) 
[24/45 |Z; 4 


(4.19) 
4,4 Surface Effectiveness. 


The total effectiveness for a finned surface may be defined as the 


ratio of the total heat transfer of the combined area of the surface and 
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fins to the heat which would be transferred if this total area were 
maintained at the base temperature T. Let "A" equal the total heat 


transfer surface including the fin area А; then 


Ho = \- PEN-A) (4.20) 


4.5 Overall Coefficient of Heat Transfer. 


Te 
ave ae p ШЕ 
Tr ME — JAVA ye A SA 
f ^ cs PES n 
Te “A 4A ZA 


Es 
Figure 7 Sketch Illustrating-One Dimensional Conduction 
and Convection through a Plane Wall 


Consider the plane wall shown above exposed to a hot fluid on one 
side and a cooler fluid on the other side. The steady state heat trans- 


fer is expressed by 
47 A. A(T«- T.)- £^ (T. — 5) - ZACS- T) (4.21) 


The heat transfer process may be represented by the resistance network 
shown and the overall heat transfer is calculated as the ratio of the 


overall temperature difference to the sum of the thermal resistances. 


$- la- — (4.22) 
А + A 


The overall heat transfer coefficient "U", is defined by the re- 


lation 


4 = (Л АД ТГоу= кам (4.23) 
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where À is some suitable area for heat flow. Then 


JE ЕЕ 58 


Now consider a plate of thickness "a" with extended surfaces or fin 


areas on either side. The overall heat transfer coefficient becomes 


| 
[accus (4.24) 
р. о 1 7 
Тон Ky (Aw/Au ) 4 (Ac/Au ) Aoc Z. 


where U is based on a unit of hot side total area including fin or ex- 
tended area. A denotes the average wall area. In the computer program 
the very small wall resistance is neglected in calculations so that 


| 
|J = 


| | 
— = lg V. 
oa Z u (Ac/A a) Moe £c 


4.6 Exchanger Heat Transfer Effectiveness 
Tu ан ТЗ cq n _ _ 


а Б (=e), 
T Per cu 


Те,дит 





Ты ОТ 


FLOW БЕЛЕТИ Teu 


Figure 8 Temperature Profiles in Counterflow 
Heat Exchanger 


The actual heat transfer in a heat exchanger may be computed by 
calculating either the energy lost by the hot fluid or the energy gained 


by the cold fluid. Referring to Figure 8, 


+" (m (e), (TÓ = h i) = (Gap (Te oor - Те, м) (4.25) 
Heat exchanger effectiveness is defined as 


actual heat transfer 
maximum possible heat transfer 
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To determine the maximum possible heat transfer for the exchanger it is 
necessary to recognize that this maximum value could be attained if one 
of the fluids were to undergo a temperature change equal to the maximum 
temperature difference present in the exchanger, ( Тим — Te, им ye 
The fluid which might undergo this maximum temperature difference is the 
one having the minimum value of /"w-Cp since the energy balance re- 
quires that the energy receíved by one fluid is equal to that given up 
by the other fluid. If the fluid with the larger value of Mm Ey goes 
through the maxímum temperature difference, the other fluid would under- 
go a temperature difference greater than maximum. The maximum heat 


transfer is therefore expressed as 


Fmax = (m ШТ qm (4.26) 


The "minimum" fluid may be either the hot or cold fluid depending on the 


mass flow rates and specific heats. 


"am (nae le ШО шыл) = (nin Cp )e (Те рот и) (4.27) 


(44. Cp 7 (Тним ~ Т, 1A ) C (Se == Те, ж) 


An expression for effectiveness in counterflow will now be derived. Re- 
ferring to Figure 8, the heat transferred through an element of area dA 


may be written 


dg = (= Ce] Та = (4). АТ 5 Cue Ты, (Ат; (4.28) 


The heat transfer could also be expressed as 


44 = (Т Был (4.29) 


Equations (4.28) and (4.29) combine to give 
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Ps Ta d UAA 
AAMA | E am (4.30) 
Tu- Te | G | Са 


Integration yields 





r í am". T Te, ES - | A Cu U A (4.31) 
Givin. Те, 7 @ e С Ц 
Or 
Cu а 
(Taout — 15, ) , pe Cc | Cu (8.32) 


l ss Te ovr) 
If Cue оя CE Cu т АТ” 








= (4.33) 
2 В QU. Cc. Ше 2, 
Tu, = Te, wad LL uc Te, od 
Then 
Та аот = Taw БТР е e” — г ғ) (4.34) 
апа 
E (n (4.35) 
Те вит E ЕУ + Є e. (Ты,ғ-Там) 
From (4.32) 
ОА 
77 = = (Ти, А) Т) — №, УУУ, s 5 № a Cw 
оа | 722 
д n7 Te, /М -е--- = (T, Are Te, 4) 
Er = s 
a \ — — 
S ANNE = = uc (> (4.36) 
l— € 8 ¿ 
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Solving for с, 


Ри еке (4.37) 


раба АШ аф 


If ic or “ - CC , the same procedure yields 


Bu e] 





= (4.38) 


e to 
[, Ce} | UA. 
ел e 


Cu 


||- Goud) UA 
== ls C m AX СРИ 


In either case 





C miM \- см А. 


= OS e Cmax | C nig 


z ам МТУ 


-— C MAX 


e 


= 
w —^ 


| Сим шығасы 
~ Cmax C SIS. 


where 


А 
ATU = ya- 


is defined as the number of heat transfer units. In general it is 


possible to express 





С = Ø (NTU, саш j Flow arrangement) 
MAX 
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vo df _. (4.39) 


(4.40) 


Ds Exchanger Flow-friction Theory. 


"uL ААА ZZ | 





Figure 9 Heat Exchanger Core Model for Pressure-Drop 
Analysis 


5.1 Core Friction. 

Friction factor is defined on the basis of an equivalent shear 
force in the flow direction per unit of wetted area. For most of the 
surfaces considered it is a combination of viscous shear and pressure 
force but no attempt is made to separate these two effects. The pres- 
sure force (form drag) is most significant in interrupted or wavy fin 
surfaces where fluid separation from these surfaces occurs, Thus 

ъ- Seu (5.1) 
or 
£ = 2-2. (5.2) 
БЕ 
It follows directly then that 


Аса 
P (5.3) 


AP- Au ок 


+q 


where < equals the average density in the core and "G" is based 
on free flow area. Curves showing friction factor as a function of 
Reynolds number are published in Kays and London (1). They are plotted 


for the same surfaces as the Colburn j-factor mentioned in the previous 
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section. 


5.2 Flow Acceleration. 
A term accounting for flow acceleration must also be included be- 
cause of the change in density as the fluid flows through the core. A 


force balance using change in momentum yields 
| = m (U L- Ua) (5.4) 


or 





mo. m(G_G\. am 
A\ 45(% = 3 Pa (4 | (988) 


The core friction term and the flow acceleration term may now be com- 


bined to give the pressure drop between points "a" and "b". 


Aue ER Ec E be en (5-6) 


5.3 Entrance and Exit Effects. 

À typical gas turbine regenerator installation usually involves a 
flow contraction at the heat exchanger core entrance and a flow expan- 
sion at the core exit. They introduce additional flow stream pressure 
drops. This problem has been investigated both analytically and experi- 
mentally by Kays (5). 

The entrance pressure drop is made up of two parts. The first, 
(Pi - Pu is the pressure drop which would occur due to a flow area 
change alone. The second, (Pi - Po? is the pressure loss due to the 
irreversible free expansion that always follows the abrupt contraction. 
First considering the pressure drop due to area change alone, 

й + ам +4 B Ua (5.7) 
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The pressure change from "Il'" to 'a" is very small relative to the total 


pressure so that 2 E - Therefore 
As. Z = A. Жа (5.8) 


Substitution of equation (5.8) into (5.7) yields 
(&-),- dg wt (i-es Qo) 
iR) - zs AW (im Ie (5.9) 
The analysis, as presented in Kays, of the pressure drop due to 
free expansion is somewhat more complex and only a brief description of 


the equations utilized in the computer program will be given here. The 


contraction loss coefficient, K. is defined as follows: 


(Б fah a % гй ктер (5.10) 


K. can be evaluated from 


= 
Cis, 


psa CSS (pé -1) 
on 


(5.11) 


K, is a function of Reynolds number and flow geometry. For laminar 


d 

flow 
К, - 1.333 (circular tubes) (5.12) 
К, - 1.390 (square tubes) 
K, = 1.430 (triangular tubes) 


For turbulent flow in circular tubes 


KA=1.09008 (45) 4.0L 694/45 +| (5.13) 

where 
—0.2 
$= 0.044 Re (5.14) 

and 

Kdisqvacey—| _ _ КА (татАМв Еу 1 _ (5.15) 

mo a a 
KaAccireucar) — | 2 KA (CIRCULAR) — | 
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Са is a function of the contraction ratio Q` and is presented graphical- 
ly for circular tubes and a one-dimensional contraction. Kays recommends 
that the one-dimensional contraction be used in evaluating triangular 
and square tubes. The curves of C. vs. G were fit with straight line 
approximations for use in the computer program and are shown in Ap- 
pendix B. 

The exit pressure rise is similarly made up of two parts. The 
first, (P, - Poy is the pressure rise which would occur due to area 
change alone. The second, (P, - Р,),, is the pressure loss associated 
with the irreversible free expansion and momentum changes following an 
abrupt expansion. Тһе first term follows from the same analysis as the 


corresponding term in the entrance pressure drop. 





š ay = же A gt n) = Ez (r) (5.16) 
The free expansion term is given by 
N a 
p, LR). 5% s = == (5.17) 
where 
Ke= 1-2 Kd oe әт (5.18) 


and Ка is evaluated using (5.12), (5.13) and (5.15). 

K. and K, are presented graphically in Kays and London for the 
circular, triangular and square tubes; however, the preceding equations 
were used for evaluation in the computer program. 


If the core friction, flow acceleration, and the entrance and exit 


effects are now combined 
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- £i л Your +. 
2% (1- otk.) (9 —\|+&А 22 (ска | (5.19) 
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a Computer Program Development. 
6.1 General Approach. 

As discussed in the introduction, the program simulates a gas tur- 
bine arrangement running under constant load conditions. The nature of 
the problem necessitates that an iterative approach be taken in solving 
for most of the desired outputs. The program is flexible in that six 
different sets of output data are obtainable depending on what type of 
information the user desires. 

Program inputs can be divided into four categories as follows: 

a. known information describing the cycle as a whole 

b. information describing regenerator core characteristics 

c. assumed quantities necessary to start the solution 

d. special inputs depending on the type of output desired 
The first three categories are referred to hereafter as general inputs. 

Program outputs can be divided into two categories as follows: 

a. special outputs depending on which set of output data is 
desired 
b. general outputs common to all sets of output data 

The program is executed largely through the use of subroutines so 
that characteristics of various gas turbine arrangements are easily 
interchangeable. All general outputs are solved for in one rather 
long subroutine called XMAIN.  Subroutine XMAIN utilizes all the other 
subroutines as the solution progresses. Subroutine XMAIN is described 
in detail later in this section; however, the general approach for the 
solution of the various unknowns in the gas turbine cycle is as fol- 
lows: 

a. Referring to Figure 1l, the temperature and pressure at 


point "2" are solved for knowing the temperature and pressure at point 
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ті", the compression ratio and the compressor efficiency. 

b. The pressures at points '5'" and "3" are computed using 
assumed values for pressure drops through the regenerator. 

c. An estimate of the temperature at point '5' is computed 
using ШЕ and the compression ratio. 

d. The temperatures at points "3" and '6'' are computed using 
the estimated temperature at point ''5'' and an assumed value for re- 
generator effectiveness. 

e. Better values are computed for the regenerator pressure 
drops using theory from section 5 and an assumed mass flow rate. 

f. Better values are computed for the temperatures at points 
"5" nana ''6'', 

g. A new value is computed for regenerator effectiveness 
using the theory from section 4 and fed back into step "4". 

h. The process is repeated until all temperatures and pres- 
sures converge. 

i. A new value of mass flow rate is now computed by comparing 
the calculated value of power output to the desired power output. 

j. The entire process is then repeated until mass flow rate 
converges. 


k. All remaining general outputs are now computed. 


6.2 General Inputs. 


a. The following inputs describe the cycle as a whole: 


T4 compressor inlet temperature 
Р, compressor inlet pressure 

P, /Ру compressor pressure ratio 
Ne compressor efficiency 
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b. 


istics: 


Җеом» 


4A P 


mm = саана: 


з 
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Te 


ЧА 
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(2 
4 


Е». 
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Үме 
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А. 


А |, 


а 
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turbine inlet temperature 

hot side regenerator exhaust pressure 
turbine efficiency 

combustor efficiency 

relative pressure drop through combustor 
net horsepower output 

fuel temperature at combustor inlet 
lower heating value of fuel 


gas constant for air 


The following inputs describe the regenerator core character- 


density of regenerator material 


cold side heat transfer area/volume between 
cold side plates 


hot side heat transfer area/volume between 
hot side plates 


cold side plate spacing 
hot side plate spacing 
cold side hydraulic radius 
hot side hydraulic radius 
cold side fin thickness 
hot side fin thickness 


cold side fin or extended area/total cold side 
transfer area 


hot side fin or extended area/total hot side 
transfer area 


plate thickness between cold and hot sides 


fin thermal conductivity 


The following are assumed input quantities: 


° 
m 
a 


mass flow rate of air 
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a A, P relative pressure drop through cold side of re- 


"A generator 

„Ае? relative pressure drop through hot side of re- 
Pe generator 
= regenerator effectiveness 


6.3 Subroutines Other Than Subroutine XMAIN 
a. Subroutine KEKC (SIGMA, RE, XKE, XKC) 

Input arguments are (7^ and Re ; computes Ke and Kc 

using equations from section 5. 
b. Function PRAND (T) 

Input argument is T ; computes Prandtl number of air from 
equation obtained by fitting curve to data given in Hilsenrath (6). 
See Appendix B for data. 

с. Function VIS (T) 

Input argument is T ; computes Ж of air from equation 
obtained by fitting curve to data given in Hilsenrath (6). бее Ар- 
pendix B for data. 

d. Function SPEC (A) 

Input argument is T 5 computes E of air. The equations 
used for computation in this and the next four subprograms are given in 
Howard (2). 

e. Function TEMR (B) 
Input argument is Pr ; computes T of air. Pr as used here 
refers to relative pressure, not to be confused with Prandtl number. 
f. Function ENTH (B) 
Input argument is T ; computes h of air. 
g. Function TEMH (B) 


Input argument is h; computes T of air. 
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h. Function PREP (B) 
Input argument is T; computes Pr of air. 
1. Function FUENTH (T) 

Input argument is Te; computes enthalpy of fuel from equation 
given in section 2. 

J. Function FRICC (RE) 

Input argument is regenerator cold side Re; computes cold 
side f from equation obtained by fitting curve to experimental data. 
See Appendix B for details on £ and j vs. Re data used in this and the 
following three subprograms. 

k. Function FRICH (RE) 

Input argument is regenerator hot side Re; computes hot side 

f from equation obtained by fitting curve to experimental data. 
1. Function STPRC (RE) 

Input argument is regenerator cold side Re; computes cold side 

j from equation obtained by fitting curve to experimental data. 
m. Function STPRH (RE) 

Input argument is regenerator hot side Re; computes hot side 

j from equation obtained by fitting curve to experimental data. 
n. Function GENFX (TU, RATIO) 
Input arguments are Сөніп Салах and NIU; computes € from 
equation developed in section 4. 
о. Subroutine XWRITE 
This subroutine contains all general output print statements. 
р. Subroutine XOPEN 


Computes open cycle outputs 
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6.4 General Outputs. 


a 


= 
V $a. 


SFC 


Wt. 


с 


һ 
NTU 


mass flow rate of air 

regenerator effectiveness 

fuel to air ratio 

specific fuel consumption 

regenerator core volume 

regenerator core weight 

compressor discharge pressure 
compressor discharge temperature 

cold side regenerator outlet pressure 
cold side regenerator outlet temperature 
turbine inlet pressure 

turbine exhaust pressure 

turbine exhaust temperature 

regenerator hot side outlet temperature 
regenerator cold side Reynolds number 
regenerator hot side Reynolds number 


number of exchanger heat transfer units 


Cmin/Cmax capacity-rate ratio 


үз 


ap 


A mu» AP 


) PL қ Pe 


6.5 Detailed Developmentof Subroutine XMAIN 


Assumptions made in the development of subroutine XMAIN are: 


fluid properties of gases leaving the combustor. 





2 
(А M ‚ сом , and 50 remain constant dur- 


S 


ing the iteration process. 
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The first step in the subroutíne is to calculate T, from T,, P 





2 1 1 
апа tas ç 
Р 

ha OT) (6.1) 
Kea d (T) (6.2) 
(ед, = еа | > (625) 
«M / ( El (6.4) 
his= Ø (Tas) (6.5) 
hz * (hihi) + hi (6.6) 
Т.- Ф(А.) (6.7) 

Next, h, and 7 re127* computed from Tye 
а = 2 (Т) (6.8) 
Пе = D (Ta) (6.9) 
Р. апа Р, are now computed using the assumed values for AT and 

2. 
c Ac P | 
G Ж 
ia E r =) (6.10) 
B-R-(aet\e (ӨСТІП) 
Р. 


R= R/E) = 


It next becomes necessary to get an estimate of T This is done by using 


5° 


a cycle with no regenerator and no pressure drop through the combustor. 
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is т 


—r 


Р Б. 


б; 


Р-Р) < 
Tes = Ø (Rees) e 

hes d Tes) (6. 

hs ha He(ha- hss) (6. 

Te Bhs) e 


T4 and T 


and assumed value of Е 


5 


Т = €(Te-h) +Th (6. 


(6. 


6 are now computed using the value just computed for T, and 


13) 


14) 


15) 


16) 


17) 


18) 


19) 


lk = Te -é (Te- Ta) o 


Temperatures and pressures at points '2", "3", "5" and "6" are now 


known or have been estimated. Better values are computed for the pres- 


sure drops through the regenerator using the theory from section 5, an 


assumed value of m and fluid properties computed from the temperatures 


and pressures at the above four points. For the regenerator cold side: 


= (Уе) (6.21) 
V= R 1А (6.22) 


D Ба (6.23) 
Ул, c= Var Us + Vs (6.24) 
2. 
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= ш T (6.25) 
А 
а E (6.26) 
: m s 
д оо (6.27) 
ОЕ о (6.28) 
Åc c= cz X Àsv (6.29) 
e = АН (6.30) 
Е 9 T...) (6.31) 
_ A4 С 
Кас — (6.32) 
on Е d (sc) (6.33) 
Ge | у А. И 
lur dv (-е tKa) H(I) t & BE E 
ба (92° —| + Ke,c | 
(6.34) 


The pressure drop through the regenerator hot side is computed in the 


same manner except that the fuel to air ratio must be included in cal- 


culating СЕ 


СЕ жаа (14- Сұ) (6.35) 


С 


where 


> (а-а (6.36) 


(К+ Kom 4 ka) 
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P, and Р, are now recalculated using the results just obtained. 


3 
= (&-R) (6.37) 
B -Qa (&-R) (6.38) 
is computed as follows: 
& = R- (SPR (6.39) 
3 Pa © ° 
TAE Rer, ( = (6.40) 
Tes = d (Pas) (6.41) 
hss = Ф (Те) (6.42) 
M la - Te (ha hes (6.43) 
Те = Ф (Ағ) (6.44) 


À new value of T; 


Using this value of T, and the assumed C  , new values for T, and T 


5 3 6 
are computed. Due. Tah, Ate Е ha; `+. ‚бү, Re. and Re, 


are then recomputed using the above values of Т; апа Te: 


A better С can now be determined using the theory from section 


фк = 6 (Кел) (6.45) 
de= P (Rac) (6.46) 
ERT (6.47) 
ЕЛЕР (6.48) 
S з c.m 
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о, е = I-A uc (1- Қос) 
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(6. 


(6. 


.58) 


.59) 


‚ 60) 


61) 


‚ 63) 


The above value of € is now fed back into equation (6.19) and the 
entire set of calculations from (6.19) through (6.63) is repeated until 
the desired degree of convergence is obtained. Many different criteria 
for determining convergence could be used; however, numerous trials in- 
dicate that the regenerator cold side pressure drop converges at a slower 
rate than any of the other computed outputs. When the value of (P, - P3)» 
computed from equation (6.34), converges to within one percent of its 
previous value the subroutine leaves the above loop and proceeds to the 
next step. 

Thermal efficiency and the true value of mass flow rate needed to 


give the desired horsepower output at that thermal efficiency are now 


computed. 
Ae Нач) 6.68) 
Via Ga 
Then 
Hes Tee rs ma f _ (6.65) 


2544 
where mais the assumed value. The above calculated value of horsepower 
is then compared with the desired value. If the two agree within .1 
percent of the desired horsepower the subroutine ends; if they do not a 
new value of m. is calculated 
‚ _ 2744 НИ 
Her i/a 44 


where HP is the desired horsepower output. This value of m is then 


(6.66) 


fed back into the subroutine where the process of calculating pressure 
drops, effectiveness, thermal efficiency and horsepower is repeated 


until horsepower converges. Convergence is normally rapid and two or 
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at most three times through the set of calculations is sufficient. 
The subroutine then proceeds to calculate the remainder of the 


general outputs. 


ссе = 2144. (6.67) 
Mek $< 
О ce 
P X Vs (6.69) 
57 MET CANET Tu doo 
Р TOTAL EN (M Ps 


6.6 Special Inputs and Outputs. 

As mentioned before, the program is flexible in that different sets 
of output data are obtainable depending on what type of information the 
user desires. The five possible sets of output and corresponding inputs 
are described below. 

a. Output type one 
In some cases the frontal area and regenerator core 
length may be fixed at some predetermined values. The user may then 
determine gas turbine performance for that frontal area 
and length using output type one. The frontal area and length along with 
the general inputs are read into subroutine XMAIN and all general out- 
puts are computed. 
Inputs: general, length, frontal area 
Outputs: general 
b. Output type two 
It may be desirable to fix the overall thermal efficiency 


at a predetermined value. There may also be some latitude possible in 
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choosing the frontal area of the regenerator providing the length re- 
mains between specified limits. Output type two is used to determine 

a frontal area for each length under consideration that will give the 
desired overall thermal efficiency. The shortest length under considera- 
tion, the longest length under consideration, the size of the increments 
desired between the shortest and longest lengths, the desired thermal ef- 
ficiency and all general inputs are read into the computer program.  Sub- 
routine XMAIN is then called with the shortest length and an initial 
value of frontal area, determined by the computer, as arguments. Frontal 
area is then incremented with length held constant until the computed 
thermal efficiency is in agreement with the desired thermal efficiency. 
The same procedure is repeated for all other lengths under considera- 
tion. 

Inputs: general, shortest length under consideration, 
longest length under consideration, size of increments desired between 
shortest and longest lengths, desired thermal efficiency . 

Outputs: general, required frontal area for each length 
considered. 

с. Output type three. 

Output type three is similar to two in that the overall 
thermal efficiency is set at a predetermined value. In this case, how- 
ever, it is assumed that there is some latitude possible in choosing the 
length providing the frontal area remains between specified limits. А 
length must therefore be determined for each frontal area under considera- 
tion that will give the desired thermal efficiency. The same computing 
procedure is followed as in output type two except that frontal area is 
held constant and length is incremented until the desired thermal effi- 


ciency is obtained. 
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Inputs: general, smallest frontal area under considera- 
tion largest frontal area under consideration, size of increments de- 
sired between smallest frontal area and largest frontal area, desired 
thermal efficiency. 

Outputs: general, required length for each frontal area 
considered. 

d. Output type four. 

In some cases it may be desirable to find the minimum 
frontal area that will allow a predetermined thermal efficiency. In 
other words make the frontal area as small as possible without giving 
consideration to length. Output type four is used to determine this 
minimum frontal area and the length that must be used with it to obtain 
the desired thermal efficiency. 

All general inputs along with the desired thermal efficiency are 
read into the program. Subrouting XMAIN is called and an initial length, 
depending on the mass flow rate and internal core geometry, is deter- 
mined by the computer. Frontal area is then incremented until the de- 
sired thermal efficiency is obtained. A new length is then selected and 
again frontal area is incremented until the frontal area corresponding to 
the new length gives the desired thermal efficiency. This process is 
repeated and frontal areas are compared until a minimum allowable front- 
al area is obtained. 

Inputs: general, desired thermal efficiency 

Outputs: general, minimum frontal area for desired 
thermal efficiency, length to be used with that frontal area 

e. Output type five. 
The user may find it helpful to determine the maximum 


thermal efficiency possible with a fixed allowable frontal area or range 
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of frontal areas. Output type five computes the maximum thermal ef- 
ficiency and the length that must be used with the frontal area under 
consideration to obtain that thermal efficiency. The frontal area under 
consideration is read into subroutine XMAIN and length is incremented 
until thermal efficiency reaches a maximum. 

Inputs: general, smallest frontal area under considera- 
tion, largest frontal area under consideration, size of increments de- 
sired between smallest frontal area and largest frontal area. 

Outputs: general, length corresponding to each frontal 
area considered for maximum thermal efficiency. 

f. Output type six. 

Output type six is utilized in determining thermal ef- 
ficiency, mass flow rate, specific fuel consumption, fuel to air ratio, 
compressor discharge temperature and pressure and turbine exhaust tempera- 
ture of the simple open cycle. This information can be used in compari- 


son with other output types in setting a standard for regenerator per- 


formance. 
Inputs: general except regenerator characteristics 
Outputs: as given above 
g. Tabular summary of special inputs and outputs 
Output 
Type One Two Three Four Five Six 
Acc I --- --- --- --- --- 
А š --- 0 І 0 І --- 
fr,min 
А --- 0 І --- I --- 
fr,max 
L I --- --- 0 --- --- 
uiu sni : ы = 7 jo 
L --- І 0 --- --- --- 
тах 
opt --- --- --- --- 0 --- 
bra 0 I I I --- 0 


GENERAL SPECIAL 


INPUTS INPUTS 
SPECIAL 
MAIN PROGRAM OUTPUTS 
KE, KC 
| SUBROUTINE SUBROUTINE 
f, j XMAIN | XWRITE 
FLUID 
PROPERTIES 
GENERAL 


REGENERATOR OUTPUTS 
EFFECTIVENESS 


FIGURE 10. FLOW CHART OF COMPUTING PROCEDURE 
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45 Results. 

The ORENDA OT-4 gas turbine was utilized as a model in the develop- 
ment and testing of the computer program. This engine meets a require- 
ment that was generated by the U. S. Army and Navy in 1960 and is dis- 
cussed in a paper by Quan (7). The design target at 600 hp on specific 
fuel consumption was .40 lbm of fuel per hp hr, which is equivalent to 
a thermal efficiency of 34.67. The pressure ratio ranges from four to 


five. The following were additional design targets: 


compressor efficiency 86% 
turbine inlet temperature 2200°R 
turbine efficiency 89% 
combustor efficiency 98% 
regenerator effectiveness 82.57 


А few additional quantities were required to complete the general inputs 
in the computer program. Their values were selected as realistically as 


possible as follows: 


compressor inlet temperature 560°R 
compressor inlet pressure 2120 Wee 
fuel inlet temperature 560°R 

lower heating value of fuel 19,000 BTU-1bm | 
regenerator hot side exhaust 2120 IDE NE 
pressure 


Six cores were selected at random for evaluation and incorporated 
into the above gas turbine arrangement. The characteristics of the 
first five shown on the following pages are from Kays and London (1); 
the last one from Pucci et.al.(8). The friction factor and Colburn 
j-factor data utilized in obtaining functional relationships for each 


core is listed in Appendix B along with the curve fitting results. The 
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cores were made of stainless steel, aluminum or nickel; however, the 
density of stainless steel was used as the input in the program for all 
cores to obtain a common basis for weight comparison. Output types one 
through six were computed for each core. The desired thermal efficiency 
of .34 was set for output types two, three and four. Required frontal 
areas were computed to match lengths of 3.0 ft through 7.0 ft in output 
type two and required lengths were computed to match frontal areas of 

3.0 pee through 7.0 fe in output type three. All core incrementing in 
output type three was done in .1 ft increments so that lengths determined 
for (єк = .34 ате only accurate to within .1 ft. The thermal efficiency 
obtainable with these computed lengths may therefore be slightly higher 
or lower than .34. The characteristics of each core along with results 
selected from the various program output types are tabulated. 

The tabular results are followed by six sets of curves which il- 
lustrate the range obtainable in overall dimensions and weight depend- 
ing on the core selected. The set of curves shown in Figure 12 illus- 
trates the lengths that must be used over a range of frontal areas to 
obtain a thermal efficiency of .34. Figure 13 gives regenerator core 
weights as a function of frontal area to obtain a thermal efficiency of 
.34 and Figure 14 is a plot of regenerator effectiveness vs. total re- 
lative pressure drop for a thermal efficiency of .34. Figure 15 il- 
lustrates the maximum thermal efficiency obtainable as a function of 
frontal area and the optimum lengths for each frontal area, i.e., the 
lengths necessary to obtain maximum thermal efficiency are shown in 
Figure 16. Core weights as a function of frontal area for maximum 
thermal efficiency are shown in Figure 17. 

A sample of the six computer output types, generated using core l, 


is shown following the curves. Output type one is shown on page 73, 
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two on pages 74 through 76 , three on pages 77 through 79 , four on 
page 80 , five on pages 81. through 83 and six on page 84. 

As mentioned above, the six cores evaluated were selected at ran- 
dom. No systematic attempt was made in finding an optimum core from 
the standpoint of external dimension or weight considerations. The 
results serve only to illustrate the usefulness of the computer program 
in selecting the best core for a given gas turbine. A few very general 
observations regarding the geometries investigated are as follows: 

1. The desired thermal efficiency was not obtainable with the 
fine mesh cores in the small frontal area range, regardless of length. 

2, Significantly longer lengths were required to produce 
the desired thermal efficiency when the coarse mesh cores were utilized. 

3. The coarse mesh cores add considerably more weight and 
volume to the gas turbine arrangement than the fine mesh cores for a 
given thermal efficiency. 

4. Figure 12 shows that required core weight decreases or 
remains essentially constant as frontal area is increased with He, = . 34 
for all cores except core 4. It is apparent that this should be the 
case after inspection of Figure 11 which reveals that required volume 
remains essentially constant as frontal area is increased for all cores 
except core 4, 

5. Total relative pressure drop vs. regenerator effective- 
ness appears to be a linear relationship for fixed thermal efficiency 
independent of internal core geometry. 

6. Fine mesh cores are much more sensitive to changes in 
external dimensions than coarse mesh cores. Shorter length increments 
than .1 ft are necessary for good accuracy when investigating fine mesh 


cores. This is apparent from Figure 13 where the points determined for 
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cores 3 and 6 are the most seriously in error. Instructions for changing 
the length increments in the computer program are given in Appendix A. 

7]. Higher maximum thermal efficiencies are obtainable 
with the coarse mesh cores, however, tbe corresponding lengths and weights 
are much greater than with the fine mesh cores. 

8. Cores 3 and 6, having approximately the same plate spacing, 
produced essentially the same results even though core 6 had perforated 
fins and 3 did not. 

9. Core 5 produced results similar to 3 and 6 because of its 
extremely high friction factor even though the plate spacing is nearer 
to that of cores 1, 2 and 4. 

10. All cores produced a substantial increase in thermal effi- 


ciency over that obtainable with the open cycle. 
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FIGURE 11. SECTIONAL DRAWING ILLUSTRATING CORE GEOMETRIES 
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Table I Characteristics sud Program Results - Core 1 


Designation, Plain plate-fin surface 10.27 T 
Plate spacing, b = 0.544 in. 

Flow passage hydraulic diameter, 4r, = .01259 ft. 
Compactness, С = 289.93 £t2/ft 

Fin metal thickness = .010 in. 

Fin area/total area = 0.863 

Fin spacing at plates, c = 0.1946 in. 


Results at thermal efficiency = .34 


Frontal Length Regenerator Weight A 

Area (£t) Effectiveness (1b) I^ / total 
3.0 desired thermal efficiency not obtainable 
319 6.2 „793 1122 „116 
4.0 4.4 .754 910 .079 
4.5 277 ‚ 738 816 ‚062 
5.0 3. 3 .731 853 .053 
2:5 2.9 .720 825 .046 
6.0 2.7 .719 838 .042 
6.5 2:29 .716 840 038 
7.0 29 „710 832 ‚036 


Results at maximum thermal efficiency 


Frontal Thermal Optimum Regenerator Weight м. 

Атеа Efficiency Length (ft) Effectiveness (1b) P Jtotal 
3.0 . 320 5.9 8765 985 . 143 
3.5 ‚343 73 . 816 1321 . 139 
4.0 ‚359 8.5 „850 1758 „132 
4,5 5971 9,7 1075 2257 .129 
5.0 ‚382 9.7 ‚885 2508 . 114 
5.5 „391 10.2 „897 2901 ‚ 106 
6.0 . 399 11.0 . 910 3413 . 103 
6.5 . 406 1127 . 919 3932 . 100 
7.0 411 1229 ‚927 4452 . 097 
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Table II Characteristics and Program Results - Core 2 


Designation, Plain plate-fin surface 12.00 T 
Plate spacing, b = 0.250 іп. 

Flow passage hydraulic ое ШЕН ж ,009412 ft. 
Compactness, (8 = 392.7 ft?/ft 


Fin metal thickness = .006 in. 
Fin area/total area = .773 
Fin spacing at plates, c = .1667 in. 
Results at thermal efficiency = .34 
Frontal Length Regenerator Weight АБ. 
Агеа (ft) Effectiveness (1b) j> total 
3.0 desired thermal efficiency not obtainable 
255 271 ‚769 607 ‚095 
4.0 2.6 ‚749 582 ‚ 071 
4.5 22 . 7 30 554 .057 
5.0 2.0 .726 560 . 049 
b 1.5 5719 554 ‚ 044 
6.0 L: 7 „#19 571 ‚040 
6.5 1.6 ,719 582 . 037 
7.0 1:2 2719 588 035 


Results at maximum thermal efficiency 


Frontal Thermal Optimum Regenerator Weight СЭ 


Атеа Efficiency Length (ft) Effectiveness (1b) total 
3.0 2332 4.0 ‚803 672 ‚149 
39 .349 4.9 .838 960 . 152 
4.0 ‚ 365 5.8 . 866 1299 . 135 
CONS 977 6,6 .887 1663 ‚128 
5.0 . 388 ne . 902 2016 ІЛ! 
9.5 .397 7.8 . 914 2402 2116 
6.0 .404 8:3 . 924 2788 SUP 
6.5 .410 8.3 . 928 2984 ‚ 101 
7.0 415 8.4 . 933 3292 ‚096 
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Table III Characteristics and Program Results - Core 3 


Designation, Plain plate-fin surface 46.45 T 
Plate spacing, b = 0.100 in. 
Flow passage hydraulic diameter, 4r, = 0.002643 ft 


уһ ft 


Compactness, £ = 1332.45 ft2/ft 


Fin metal thickness = .002 in. 
Fin area/total area = .837 
Fin spacing at plates, c = .0431 


Results at thermal efficiency = .34 


Frontal Length Regenerator Weight AP 
Area (ft) Effectiveness (1b) 7598 
cotal 
3.0 desired thermal efficiency not obtainable 
329 desired thermal efficiency not obtainable 
4.0 0.5 . 819 206 ,139 
4,5 0.3 ‚760 139 , 084 
5.0 0.3 179 154 ‚077 
925 0.2 JL 113 .055 
6.0 222 «730 125 (OST 
6.5 0.2 ‚765 134 . 048 
7.0 0.2 2778 144 .045 


Results at maximum thermal efficiency 





Frontal Thermal Optimum Regenerator Weight ГА Б 
Area Efficiency Length (ft) Effectiveness (1b) \ р 
total 

2.0 . 303 0.5 „769 154 ‚ 186 
3.5 2423 0.6 .820 216 . 185 
4.0 ‚ 342 0.6 ‚841 247 „162 
4.5 2355 0.7 .872 326 . 164 
5.0 . 368 0.7 , 884 360 . 148 
5.5 ‚317 0.8 ‚ 905 452 .150 
6.0 ‚389 0.7 „904 432 JE | 
6.5 . 397 0.8 . 920 535 . 124 
7.0 ‚403 0.9 „933 648 ‚126 
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Table iV 


a 


a 


Designation, Plain plate-fin surface 6.2 
Plate spacing, b = ‚405 in. 


Flow passage hydraulic d 4r, = 0.0182 ft 


Compactness, “ = 204.0 ft2/ft 
Fin metal thickness = 0.010 in. 


Fin area/total area = .728 

Fin spacing at plates, c = .322 in. 
Results at thermal efficiency = .34 

Length Regenerator Weight 

(ft) Effectiveness (1b) 
6 2 .756 887 
548 ‚736 968 
5T „726 1087 
5.6 2721] 1202 
5.4 ‚716 1287 
5.2 ‚ 714 1364 
4.9 ‚712 1402 
&.6 ‚210 1426 
4,3 .708 1435 


Results at maximum thermal effíciency 


Thermal Optimum Regenerator 
Efficiency Length (ft) Effectiveness 


895 119: 9 ,848 
‚71 133 .869 
‚382 16.2 . 886 
.391 15:5 ‚897 
2380 20.7 ‚907 
‚405 21.4 .912 
411 22.4 .920 
.415 dee? ‚ 923 
.420 2 2/79 ,929 
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Characteristics and Program Results - Core 4 


D 


(As 


.082 
‚ 062 
‚Оз? 
‚046 
‚042 
. 038 


° 


Weigh 
(1b) 


1559 
2219 
3090 
3969 
4935 
5612 
6551 
7035 
7843 


035 
033 
031 


t 


) total 


. 137 
. 148 
.116 
. 106 
099 
‚091 
‚086 
‚ 079 
‚075 


Table V Characteristics and Program Results - Core 5 


Designation, Wavy-fin plate-fin surface 11.5-3/8W 
Plate spacing, b = 
Flow passage hydraulic diameter, Grą = 0.00993 ft. 
Compactness, © = 345 £r2/fr3 


Fin metal thickness 
Fin area/total area 


AS TE. 


Fin spacing at plates, 


Results 
Frontal Length 
Area (ft) 
3.0 desired 
325 desired 
4.0 desired 
mrs 0.9 
5x0 0.8 
5.5 O. 7 
6.0 T 
6.5 0.6 
7.0 0.6 
Results 
Frontal Thermal 
Area Efficiency 
B. 0 „301 
B. 9 2371 
5.0 338 
ЕЭ 892 
5.0 . 364 
B. 5 574 
6.0 „383 
6.5 „391 
7.0 .398 


= „010 1ш. 
=й 822 
соч .087 іп. 
at thermal efficiency = .34 
" En CEE 
Regenerator Weight [М 
Effectiveness (1b) {> 
total 
thermal efficiency not obtainable 
thermal efficiency not obtainable 
thermal efficiency not obtainable 
„7/72 235 .092 
396 331 29/9 
‚739 319 ‚060 
Nn 347 ‚ 054 
SEI 323 . 046 
2720 347 ‚ 043 


at maximum thermal efficiency 


Optimum 


Length (ft) Effectiveness (1b) — 


N 2 N) N) — = — oO 
мл ш) Кэ СО + C Ww ке wo 


Regenerator Weight A \ 
[2 


. 142 223 . 168 
„788 319 .161 
2921F 430 2152 
.854 596 2153 
.866 703 ‚137 
‚888 910 ‚131 
2200 1092 DOT 
2:207 1237 ‚120 
„916 1448 TIG 


*Weights are low since they were based on core length not flow 


passage length 


total 


Table VI Characteristics and Program Results - Core 6 


Designation, Perforated-fin, plate-fin surface, Solar #2 
Plate spacing, b = 0.100 in. 

Flow passage hydraulic diameter, ат, = ,0046272 ft. 
Compactness, (5 = 673. 143 ft2/f£c? 

Fin metal thickness - .0047 іп. 

Fin area/total area = .645 

Fin spacing at plates, c = .075 in. 


Results at thermal efficiency = .34 





Frontal Length Regenerator Weight AP | 
Area (ft) Effectiveness (1b) = 
total 
S desired thermal efficiency not obtainable 
329 desired thermal efficiency not obtainable 
4.0 desired thermal efficiency not obtainable 
4.5 desired thermal efficiency not obtainable 
5.0 О у ‚768 523 . 089 
D 0.6 DOS 493 .074 
6.0 0. 5 ‚ 749 449 ‚061 
525 0.4 ‚725 389 -051 
7.0 0.4 . 741 419 . 048 
Results at maximum thermal effíciency 
Frontal Thermal Optimum Regenerator Weight AP 
Area Efficiency Length (ft) Effectiveness (1b) 2 
total 
3.0 2257 0.8 . 7,22 359 .186 
3.5 2305 1240 ‚766 523 ‚179 
4.0 ‚324 1-1 ‚793 658 ‚ 164 
&.5 ‚338 12 „822 808 ‚ 155 
5.0 ‚958 1072 ‚840 897 ‚136 
2 9 ‚365 1.4 ‚871 131 . 140 
6.0 2375 1.5 ‚889 1346 . 136 
D. . 386 172 .901 1453 125 
7.0 ‚ 3993 1.6 ,913 1675 149 
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FIGURE 14. 
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FRONTAL AREA = 3.0 OPTIMUM LENGTH = 5.9 
MASS FLOW THERMAL FUEL TC AIR REGENERATOR DELTA P/P 
RATE EFFICIENCY RATIO EFFECTIVENESS VOLUME WEIGHT TOTAL 
23965. 0.320 0.0105 0.765 17.7 915. 0.143 
P2 P3 P4 P5 T2 T3 T4 T5 T6 
58.89 58.46 57.29 16.66 873. 1506. 2200. 1700. 1067. 

FUEL CON. NG. COLO NG. HOT? UNITS CMIN/CMAX 
0.418 2891. 2653. 3.09 0.966 
FRONTAL AREA = 3.5  GPTIMUM LENGTH = 7.3 
MASS FLOW THERMAL FUEL TG AIR REGENERATOR DELTA P/P 
RATE EFFICIENCY RATIO EFFECTIVENESS VOLUME WEIGHT TOTAL 
23799. 0.342 0.0099 0.816 25.5 1321. 0.139 
Р2 P3 P4 P5 T2 T3 T4 T5 T6 
58.89 58.46 57.29 16.57 873. 1547. 2200. 1698. 1024. 
FUEL CON. NO. COLO NO. HOT UNITS CMEN/CMAX 
0.392 2434. 2280. 4.19 0.972 
FRONTAL AREA - 4.0 OPTIMUM LENGTH = 8.5 
MASS FLOW THERMAL FUEL TO AIR REGENERATOR DELTA P/P 
RATE EFFICIENCY RATIO EFFECTIVENESS VOLUME WEIGHT. TOTAL 
23565. 0.359 0.0095 0.850 34.0 1758. 0.132 
P2 P3 ра P5 T2 T3 T4 T5 T6 
58.89 58.48 57.31 16.45 873. 1572. 2200. 1696. 996. 
FUEL CON. NO. COLO — NO. HOT. — UNITS CMIN/CMAX 
0.373 2094. 1989, 5.31 0.976 


c8 


FRONTAL AREA = 4.5 OPTIMUM LENGTH - 9.7 
MASS FLON THERMAL FUEL TU AIR REGENERATOR 
RATE EFFICIENCY RATIO EFFECTIVENESS VOLUME WEIGHT 
23475. 0.371 0.0092 0.875 43.6 2257. 
P2 P3 P4 P5 Т2 T3 Та T5 
58.89 58.49 57.32 16.41 873. 1592. 2200. 1695. 
PUEL CON. — NO. COLO NO. HOT — UNITS" CMIN/CMAX 
922220 1845. 1770. 6.53 0.978 
FRONTAL AREA = 5.0 UPTIMUM LENGTH = 9.7 
MASS FLOW THERMAL FUEL TC AIR REGENERATOR 
RATE EFFICIENCY RATIO EFFECTIVENESS VOLUME WEIGHT 
22891. 0.382 0.0092 0.885 48.5 2508. 
P2 P3 ра P5 T2 T3 T4 T5 
58.89 58.55 57.38 16.14 873. 1594. 2200. 1689. 
PUEL CON. — NOS COLD — NOS HOT? — QONITS смімисмах 
OS 1618. 1558. DET. 0.979 
FRONTAL AREA = 5.5 OPTIMUM LENGTH -10.2 
MASS FLOW THERMAL FUEL TO.AIR REGENERATOR 
RATE EFFICIENCY RATIO EFFECTIVENESS VOLUME WEIGHT 
2260F 0.391 0.0091 0.897 56.1 2901. 
P2 P3 P4 P5 T2 T3 T4 T5 
58.89 58.58 57.41 16.02 873. 1602. 2200. 1586. 
PUEL CON. NOs COLD NO. HOT — ONITS  CMIN/CMAX 
0.342 1450. 1403. 8.04 0.981 


DELTA P/P 
TOTAL 


0.129 
T6 
945. 


DELTA P/P 
TOTAL 


0.114 
T6 
961. 


DELTA P/P 
TOTAL 


0.106 
T6 
957. 
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APPENDIX A 
Computer Program Usage 

The computer program is written in FORTRAN IV language described in 
(9). The main program and all subroutines were written using double 
precision in order to obtain the best possible accuracy. This is parti- 
cularly essential in subroutines where least squares curve fitting is 
utilized to determine a functional relationship from a set of data points. 

Required computer storage capacity is small as no subscripted vari- 
ables or arrays were introduced while writing the program. Essentially, 
the only storage requirement is for program compilation. 

Approximately three minutes of computer time on an IBM System/360 
is required for each output type described in section six with the excep- 
tion of output types one and six which require only about 30 seconds. 

a. Subroutine XMAIN 
The main program is designed for sizing counterflow heat 
exchangers only. Subroutine XMAIN is general however and could be uti- 
lized for numerous other flow arrangements provided that function GENFX 
(TU, RATIO), called by subroutine XMAIN, is changed to a regenerator ef- 
fectiveness function corresponding to the desired flow arrangement. 
b. Subroutines Dependent on Internal Core Geometry 
The relationships of friction factor and Colburn j-factor 

vs. Reynolds number must correspond to the core geometry being investi- 
gated. These relationships are fed into the program as function sub- 
programs FRICC(RE), FRICH(RE), STPRC(RE) AND STPRH(RE). An example is 
shown on page 111. Only the cold side functions are shown as the same 
core was used for both hot and cold sides in all cases. 

Subroutine KEKC (SIGMA, RE, XKE, XKC), K, and K. as a function of 


g~ and Re, depends on the shape of the core flow passages. Subroutines 
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were written for square and triangular passages and are shown on pages 
105 and 106. 
c. Data Cards 
All general inputs for output types one through five 
are read into the computer program from the first four data cards. Ten 
columns are allotted for each input. Each of the general inputs must 
be of floating point form. The inputs with corresponding units are 


listed below in the proper order: 


compressor inlet temperature ^R 
compressor inlet pressure esis 
compressor pressure ratio 

compressor efficiency 

turbine inlet temperature °R 
turbine efficiency 

assumed regenerator effectiveness 

fuel inlet temperature °R 
combustor efficiency 

lower heating value of fuel BTU-1bm ! 
hot side regenerator exhaust pressure wt scd 
horsepower output hp 
assumed relative pressure drop 

through regenerator cold side 

assumed relative pressure drop 

through regenerator hot side 

gas constant for air ft-1b£(1bm-°R) ~! 
regenerator plate thickness in 

cold side transfer area/volume E 
between cold side plates ft 
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hot side transfer area/volume ft-lbf(Ibun-*n) 1 
between hot side plates 


cold side plate spacing in 
hot side plate spacing in 
cold side hydraulic radius ft 
hot side hydraulic radius ft 
mass flow rate lb HER 
relative pressure drop through 
combustor 
cold side fin thickness ft 
hot side fin thickness ft 
fin thermal conductance BTU- (hr-ft^-*F/£t) | 
cold side fin or extended area/ 
total area 
hot side fin or extended area/ i 
total area 
3 


density of regenerator material lbm-ft 
If output type six is desired, only the first 12 general inputs are 
necessary. The third and fourth data cards must be placed in the deck; 
however, they may be left completely blank. The data cards utilized in 
generating the results for core 2 are shown on page 95. 

The fifth data card is used to designate the type of output desired. 
The fixed point integer "l'' in column one designates output type one; 
the fixed point integer "2" in column one designates output type two 
ессе. 

The sixth and last data card depends on the type of output designat- 
ed. In all cases numbers must be in floating point format. Ten columns 
are allowed for each input beginning in column one. Units are either ft. 
or NÉ corresponding to length or area. The inputs for each type out- 


put are listed below in the proper order: 
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output type one: length 
frontal area 
output type two: shortest length 
longest length 
length increment size 
thermal efficiency 
output type three: smallest area 
largest area 
area increment size 
thermal efficiency 
output type four: thermal efficiency 
output type five: smallest area 
largest area 
area increment size 
output type six: no data card 
d. Incidental Instruction in Program Usage 
There are several potential program trouble areas that 
the user should be aware of when investigating various core geometries. 
First of all, it is important to keep in mind the range of Reynolds 
numbers over which the functional relationships of friction factor and 
Colburn-j factor have been obtained. Frontal areas must be chosen so that 
the hot and cold Reynolds numbers fall in this range. This is of course 
largely dependent on internal regenerator core geometry and mass flow 
rate. Secondly, experience has shown that the iteration process diverges 
when the total relative pressure drop through both sides of the core and 
combustor exceeds approximately 30%. This pressure drop is mainly a 


function of regenerator length and core plate spacing. In general, shorter 
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lengths are necessary as plate spacing is decreased. 

Main program statements 061, 125, 126 and 141 designate the initial 
frontal area and length used in incrementing out the desired values for 
output types two and four. These statements, as shown in the program 
printout, are set for a 600 hp gas turbine with a regenerator having a 
core plate spacing of .2 inches or greater. For cores with plate spacing 
smaller than .2 inches or exceptionally high friction factor it may be 
necessary to shorten the length in statement 125 and to increase the 
frontal area in statements 061, 126 and 141 to get the incrementing pro- 
cess started properly It may also be necessary to increase the frontal 
area in the above statements for larger gas turbine arrangements. 

Main program statements 62, 92, 125, 140, 153 and 177 designate the 
length and area increments to be used in the main program output types. 
They are presently set for .1 ft and .1 £t? but may be changed as neces- 
sary if better accuracy is desired. 


e. Listing of Program Variables 


А1 dummy variable used in first read statement 
A2 dummy variable used in first read statement 
AFINC cold side fin area/total area 
AFINH hot side fin area/total area 


AFLOWC cold side free-flow area 


AF LOWH hot side free-flow area 


AFRINC incremental frontal area 

AFRON frontal area 

AFRON1 smallest frontal area considered 
AFRON2 largest frontal area considered 
AHTC cold side heat transfer area 
AHTH hot side heat transfer area 


ALFAC 


ALFAH 


BC 


BETAC 


BETAH 


BH 


CAPC 


CAPH 


CEFF 


COBJC 


COB JH 


COMEFF 


CPC 


CPH 


CRATIO 


DEN 


FC 


FCOND 


FH 


FINEFC 


FINEFH 


FINTHC 


FINTHH 


GC 


GENEFF 


GH 


H1 


total transfer area of cold side of 
regenerator/total regenerator volume 


total transfer area of hot side of 
regenerator/total regenerator volume 


cold side plate spacing 

cold side transfer area/cold side volume 
hot side transfer area/hot side volume 
hot side plate spacing 

cold side capacity rate 

hot side capacity rate 

compressor efficiency 

cold side j factor 

hot side j factor 

combustor efficiency 

specific heat of cold side air 
specific heat of hot side air 
compression ratio 

density of core material 

cold side friction factor 

fin thermal conductance 

hot side friction factor 

cold side fin efficiency 

hot side fin efficiency 

cold side fin thickness 

hot side fin thickness 

cold side mass flow rate per unit area 
regenerator effectiveness 

hot side mass flow rate per unit area 


enthalpy of air at compressor inlet 


ИТ 


H2 

H3 

H4 

H5 
Н25 
Н55 
НЕАТС 
НЕАТН 
HETVAL 
HF 
HP 

Pl 

r 

P3 

P4 

P5 

P6 
P23 
P56 


P23R 


P34R 


P56R 


PR1 


PR2S 


PR5S 


PRA 


PRC 


compressor discharge enthalpy 


combustor inlet enthalpy 


turbine inlet enthalpy 


turbine exhaust enthalpy 


compressor discharge ideal enthalpy 


turbine exhaust ideal enthalpy 


cold side heat transfer coefficient 


hot side heat transfer coefficient 


lower heating value of fuel 


fuel enthalpy 


horsepower output 


compressor inlet pressure 


compressor discharge pressure 


combustor inlet pressure 


turbine inlet pressure 


turbine exhaust pressure 


regenerator hot side exhaust pressure 


pressure drop through regenerator cold síde 


pressure drop through regenerator hot side 


relative pressure drop through regenerator 


cold side 


relative 


relative 
hot side 


relative 


relative 


relative 


relative 


pressure 


pressure 


pressure 


pressure 


pressure 


pressure 


cold side Prandtl 
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drop through combustor 


drop through regenerator 


- 


at compressor inlet 
at point 2s. 

at point 5s. 

at turbine inlet 


number 


PRH 


RADC 


RADH 


RATIO 


REC 


REH 


RELP 


КЕА 


SFC 


SIGMAC 


SIGMAH 


STC 


STH 


SUREFC 


SUREFH 


SVOL 


Т1 


T2 


T3 


T4 


T5 


T6 


T2S 


T5S 


TCM 


TEFF 


hot side Prandtl number 

gas constant for air 

cold side hydraulic radius 
hot side hydraulic radius 
capacity rate ratio 

cold side Reynolds number 
hot side Reynolds number 
total relative pressure drop 
fuel to air ratio 

specific fuel consumption 


cold side free flow area/cold side 
frontal area 


hot side free flow area/hot side 
frontal area 


cold side Stanton number 

hot side Stanton number 

cold side surface efficiency 

hot side surface efficiency 

solid core volume 

compressor inlet temperature 

compressor discharge temperature 
combustor inlet temperature 

turbine inlet temperature 

turbine exhaust temperature 

regenerator hot side exhaust temperature 
temperature at point T 


2s 


temperature at point Tss 
mean cold side temperature 


turbine efficiency 
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TF 


THEFF 


TU 


V2 


УЗ 


V5 


V6 


VCM 


VIC 


VIH 


VOL 


XKCC 


XKCH 


XKEC 


XKEH 


XL 


XL1 


XL2 


XLINC 


XMDOT 


fuel temperature 

thermal efficiency 

mean hot side temperature 

number of heat transfer units 
overall heat transfer coefficient 
cold side inlet specific volume 
cold side outlet specific volume 
hot side inlet specific volume 
hot side outlet specific volume 
mean cold side specific volume 
mean hot side specific volume 
cold side viscosity 

hot side С 

regenerator volume 

regenerator plate thickness 
regenerator weight 

cold side contraction coefficient 
hot side contraction coefficient 
cold side expansion coefficient 
hot side expansion coefficient 
regenerator length š 
shortest length considered 
longest length considered 

length increment 


mass flow rate 
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APPENDIX B 
Curve Fitting 
It is frequently necessary to represent by a functional relation- 
ship data that is given in the form of a set of x-y poínts. One may 


want to know the value of "y" corresponding to some "x" that lies between 


two tabulated "x" values or may want to find the "y" value corresponding 


to an "x" value outside the range of the tabulated "x" values. The 
problem then is to find a curve that approximates the given data with 
sufficient accuracy. 

Before the computer program could be fully developed, various func- 
tional relationships of the type mentioned above had to be obtained. 
Many of the necessary relationships were available in literature; how- 
ever, seven of them were not and had to be obtained by some appropriate 
curve fitting method. The "least squares" method was utilized in fitting 
six of the curves and straight line approximations were used on the re- 
maining one. It was decided to use the least squares method because of 
its simplicity and the smoothness of fit that is realizable in most 
cases. 

a. Least Squares Theory. 
The deviation at a data point is defined as the difference 

between the tabulated y-value and the y-value computed from a functional 


relationship. The least squares method asks for a minimum of 


N ; 4 
Ж Ur) (8-1) 
CsI | 


where N equals the number of data points, y' denotes a value computed 
from the functional relationship and y denotes the actual tabulated value. 
This leads to equations that in many caseg of practical interest are 


linear and, in principle at least easy to solve, If it is assumed that 
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the functional relationship desired is of the form 
K 
4o aera хаха + $ с p AX (8-2) 


am aj, а, каска, а, musk pe decermined to minimize 


S= (età ge] 


Taking the derivative of S with respect to а, 


(B-3) 


and setting it equal to zero 
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This gives one equation. Бу taking the derivatives of S with respect to 


through а,» another "k'" equations are 


Og = x 42 
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The k+l equations can then be solved for the k+l coefficients. 

The value selected for k can be any value less than N. For best 
results, however, k should be made only large enough to give desired 
accuracy. As k is increased the computer round-off errors inherent in 
solving the set of simultaneous equations becomes more significant and 
accuracy is lost not gained. The generated curve may also tend to fluc- 
tuate between data points if k is too large, making a smooth fit impos- 
sible. 

b. Gauss Elimination. 
The method of Gauss elimination was used for solving the 
simultaneous equations mentioned above. A system of three equations in 
three unknowns will be used to illustrate the procedure because the 


nomenclature becomes somewhat complicated for the general case оЁ "п" 


equations and "п" unknowns. Let ху, Х and x, be the unknowns and the 


2 3 


equations are 


Ан Х, + а. Х +A X = £ 
(B-6) 


А a 


Ars Xi +@й,4Х+@з,Х;= Lz 


The equations must be arranged so that 811 is not zero. Next define a 


multiplier, 


я 
y 


| 


M= (B-7) 


© 


The first equation is then multiplied by m, , and subtracted from the 


second equation. The result is 
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(Avy may) Xi + (Arr mie Ka od (Ars = m, Aia) Хз = an, Ay (B-8) 


But 
Ar)" Mri n ire 4, =O 
so ху has been eliminated from the second equation. 
Define: 
с = Avi Mi Air (B-9) 
(es = Дуу м. аз (B-10) 
— ~ mt, (8-11) 


then equation (B-8) becomes 


i / 
Аха, ee (B-12) 


If a multiplier is now defined for the third equation, 
А31 

fn m л В-13 

15 a ( ) 


and the first equation is subtracted from the third, the result is 


/ / 
Аҙ. X, X433 X3 — ә (B-14) 
where 

/ 
A us = 2117 ЖА; Air (B-15) 
( 33! - 433 — 13 C. 13 (B-16) 

/ 
Du 45 — ж; е, (В-17) 
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The resulting three equations in three unknowns are 


Au X, T OX, + 213 Х; = 4, 


/ = 
Qi X +A, X; = f. (5-18) 
/ / / 
Ai X + Ду Хз = 4; 


If x, is now eliminated from the last equation using the same procedure 


2 
the following is obtained: 


Qu X, даъ + AnX = 4, | 
ans Xi T da Xs 7 fa^ (B-19) 
аз = 


It is now a straight forward process to solve for x to substitute that 


3? 
result into the second equation to find хо, and finally to substitute in- 
to the first equation to find x): 


/ 





%- A // 
22 


/ / 

EO AE ES (8-20) 
Arr 

Ха (Са ~fi Xi Д1 Хэ 


Д | 


c, Curve Fitting Computer Program 
The computer program shown on the following pages was 
utilized in fitting the necessary curves. The main program manipulates 
the data point information into the equation form shown in part "a" of 
this appendix. Subroutine GAUSS 2 is then called and the equations are 
solved simultaneously using the method described above (10). The coeffi- 


l E а, determined for a polynomial of the form y = a, + 


1 1 
ax + a x° о hut are printed out along with a column of x-y 
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points generated from the polynomial. These points can then be compared 
with the original data points to determine whether or not the fit is 
sufficiently accurate. 

À maximum of fifty data points can be read into the computer. For 
best results, they should be as evenly spaced along the "x" axis as pos- 
sible. А special problem arises when using the program to fit Colburn-j 
and friction factor vs. Reynolds number data. This data is normally 
plotted on logarithmic coordinates and is therefore not evenly spaced on 
linear coordinates. To remedy this situation, the natural logarithms of 
the "x" and "у" values can be taken before starting the curve fitting 
procedure. The resulting polynomial is of the form Іп(у) 21 ar a, (1n x) + 
a, (1n Di .......л(іп ju Anti-logarithms can then be taken to deter- 
mine the desired "y" value. These additional steps are incorporated into 
the computer program as it is shown. If the user desires to fit data 
plotted on linear coordinates, cards 150, 160, 510, and 520 must be re- 
moved and card 530 must be changed to read: 

WRITE (6,55) X1(IQ), Y1(10) 

Coefficients shown in parts 1 and 2 of Table VII are for polynomials 
of the first form shown above. The remaining coefficients in Table VII 
are for polynomials of the form just described. 

When utilizing the program, the following data cards are necessary: 

(1) The number of data points to be used in fitting the 
curve is punched into the first three columns on the first data card. 
The order of the polynomial desired is punched into the second three 
columns of the first data card. Both of these numbers must be in fixed 
point integer form. The order of the polynomial can be any number up to 


one less than the number of data points. 
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(2) The x-y values of the data points are punched into 
the next N data cards, one data point per card. The x-value is punched 
into the first ten columns; the y-value is punched into the second ten 
columns. Both numbers must be of floating point form. 

(3) As mentioned above, the printout includes a column 
of x-y points generated from the computed polynomial. The first point 
is computed using the smallest x-value and the last point is computed 
using the largest x-value of the original data points. The size of the 
desired increments between these two data points is punched into the next 
data card in the first ten columns. Floating point form must be used. 

(4) If the user desires to fit more than one set of 
points he can place the next set of data cards, for those points, be- 
hind the set just described. 

(5) The next and last data card must contain the integer 
2 punched into column three. This signals that all data point sets have 
been fit and stops the program. 

d. Curve Fitting Results 

A sample of least squares curve fitting is shown in Figure 
18. All data points used in curve fitting and the resulting polynomial 
coefficients are shown in Table VII. 

The curve fit by straight line approximations is shown in Figure 

19, with the straight lines drawn in. Two straight lines were used, one 
covering the range,.3< (Q^ < .55 the other covering the range 
.55< O-<.75. This approximation was used rather than the least squares 
theory because small errors in C. cause virtually no error in the comput- 


er program results. The curve is from Kays (5). 
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Table VII Data Points Used in Curve Fitting and Resulting 


Polynomial Coefficients Used in Prograia Subroutines 


l. Dynamic Viscosity of Air ac Atmospheric Pressure from 
Hilsenrath 16). У = 1716 x poe 


94 
ТУК fW 


— me 


POolynomial Coefticiencs 


180 ORO as А) = + 4132582012 D-01* 
550 ТОА A? me 27355893209 ПЕ02 

900 113555 A3 -.1318?733093 D-OS 

1260 172406 $4 „2538864121 D-12 

1620 | 85: 12215554] 9857 ee 

1980 205262 Аб = =, 27152187450 0-15 
2240 22974 A? 2117995215 D-18 
2700 3.065 AB -.3862860359 D-22 

2060 3.290 АЗ 2782105119 D-26 
seo Beo 419: 61632858289 0-32 


и. Prandtl Number of Air fromen ea en (6). 





T^R Prandtl No. rolynomial Coefficients 
180 0,7709 Al = .7865630968 9-00 
350 Ú 729 42 = -.2871582338 0-04 
540 0.708 АЗ = - 3666074673 9-06 
? 20 0.689 А+ = -, 1373281158 0-09 
900 0.680 А5 = .1569904076 Р-11 
1080 0.680 А6 = -.1973361281 р-14 
1260 0.684 А? .1021204710 В-17 
1440 0.689 A8 -.204725292à D-21 
1620 0.696 Ад « -,1538501631 D 25 
1800 2702 А1О = „7960530427 D-29 


3. Core 1 (Plain plate-Ei® surface 102797) Friction Paetor 
fron Xays and London (1 


Reynolds Friction 
Number = Factor tolynomial Coetficients 
360 „066 А} = ,10624542201 D«03 
400 S A2 - -,/242808930 D402 
500 „0542 AS = 1760821881 PTO2 
700 „032 Аб ж -,1171044519 р491 
900 ‚026 А5 = -,2766351354 D+00 
1250 ,0185 Аб = 6028901137 р-01 
1505 ‚ 9165 АТ = -, 1810176966 р-02 
2000 20135 AB * -. 2728780735 D-03 
3000 КОШ д9 = ,1143227913 р-03 
4000 . 0095 AiO - -.6409201873 D-05 
6000 20092 411 = .1272474507 р-06 


*А] = -.4132588912 х 10-1 


4. 


Reyaolds 
Number 


300 
400 
500 
700 
900 
1250 
1500 
2000 
3000 
4000 
6000 


35 


Reynolds 


Number 


200 
300 
500 
700 
1000 
1500 
2000 
3000 
4000 
6000 
3000 


б. 


Reynolds 
Number 


200 
300 
500 
700 
1000 
1500 
2000 
3000 
%000 
6000 
8000 


Core 1 (Plain plate-fin surface 10.27 T) Coburn-j 
Factor from Kays and London (1) 


Coburn- j 
Factor 


‚016 
20/153 
zl» 
. 0087 
20072 
„006 
2905175 
‚0045 
‚ 00375 
‚00345 
‚0032 


Polynomial Coefficients 


Al 
A2 
A3 
AG 


- 9478036457 
25195775000 
.1236662988 
. 2136908795 
- 1534935252 
22255005559 
. 3614040571 
291505173943 
41222790555 
. 1473493662 
.9202049786 


D+0 1 
D+01 
D+01 
D+00 
D-r00 
n-01 
D-03 
D-03 
D-04 
D-05 
D-07 


Core 2 (Plain plate-fin surface 12.00 T) Friction 
Factor from Kays and London (1) 


Friction 


Factor 


‚080 
8057 
1057 
1927 
2919 
15 
2013 
0115 
. 0112 
"0095 
‚0085 


Polynonial Coefficients 


А1 
А2 
АЗ 
AS 
A5 
А6 
А? 
А8 
А9 
А10 
А11 


uH ow H H wu Wd ot 


B H H N 


229556554237 
4366648607 
. 7695241634 
.3558930941 
.5846569710 
‚2516782205 
.3534066422 
.37196229960 
252119225359 
.1641370472 
. 5085864719 


Core 2 (Plain plate-fin surface 12.00 T) 
from Kays and London (1) 


Coburn-j 


Factor 


.016 
012 
‚008 

. 0064 
‚0049 
.0038 
200555 
.0036 
70035 
200325 
. 003 


D-+03 
D+02 
D+02 
D+02 
D+01 
D-+00 
D-01 
D-02 
D-04 
D-04 
D-06 


Coburn-j Factor 


Polynomial Coefficients 


Al 


l H m Hu H N 


ШІСІ IH 


. 1678025363 
. 1147089449 
. 2302525106 
99307793908 
‚8790472874 
"7256951992 
.5905850854 
. 6125863884 
.9700663367 
. 1394667289 
.4807945572 
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D--00 
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7. Core 3 (Plain plate-fin surface 86,45 T) Friction Factor 
from Kays and London 1: 
Reynolds Friction 
_Number Factor Polynomiai Coefficients 
300 .059 Al = .5218554745 рн2 
400 .045 2 = 3353326320 3405102 
500 no АЗ = „8259113372 95] 
700 10265 А& = -.180226?981 D+O1 
900 2020 45 = .,4388704605 D-400 
1250 20153 Аб = -,4819376576 В-01 
1500 0152 А) = -,3956062581 D-02 
2000 0120 А8 = .7437300503 0-03 
3000 0115 А9 = .3659712646 р-04 
4000 OZ ALG = -.1134921023 D-04 
$009 ОО Ali = .5000930552 0-05 


8, Gore 3 {flain plate~fin surface 46.45 Т) Coburn-j Factor 
from Kays and London (1; 


Reynolds Coburm- j 
Number | grec tor Polynomiai Coefficients 
360 ,0135 Al = -.9179528233 D+O1 
829 „ОТОУ А2 = „2559082428 0+01 
500 . 0086 АЗ = -,1138319133 2501 
708 ‚ 0055 AS & .131231/5 5M DO] 
960 -0053 45 = -.5200720093 D+00 
1250 ,0041 AG = ,6767637499 D-01 
1560 ‚ 0036 АЈ = .2230241152 0-02 
2000 . 60295 А8 = -.1010311829 0-02 
3060 ,0026 A9  .1205538219 D-04 
4000 20025 410 = .8449625214 D-05 
6000 100225 All = -.5529482542 1-06 


9. Core ^ (Piain plate-fi^ surface 6.2) Friction Factor 
from Kays and London 11 ! 


Reynolds Friction 
Number. Factor. Polynomial Coerficients 
600 ‚027 81 = .4512980279 р+05 
300 ОРУ A2 * -.3083071127 р405 
1250 , OS АЗ = 77270459904 D+04 
2600 . 010 Аб ғ -.7810141046 0503 
2300 ‚0098 45 = -.6317179578 D+O 1 
2700 20093 Аб ж .3880647509 р+01 
3000 . 0092 27 = .6306627426 2-01 
4000 ‚ 00875 А8 = ‚3219487556 0-01 
5000 ,0081 А9 = -,1290508496 D-O1 
8000 . 00775 ALO = ‚1130983538 р-02 
10000 20099 All - -,3189159559 р-04 
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10. Core 4 (Plain plate-fin surface 6.2) Coburn-j Factor 
from Kays and London (l. 
Reynolds Coburu- j 
Number Factor Polynomial Coefficients 
600 20072 Al = .7203819966 D+05 
800 ‚0059 A2 = -.4896069217 D+05 
1250 20059 АЗ = .1220436223 D+05 
2000 ‚ 0032 A4 - -.1162973846 D-404 
2300 . 003 А5 = -.1005649404 D+02 
2700 .0031 Аб = .6019720585 0501 
3000 20035 A7 = .9239757582 9-01 
4000 . 0034 А8 - .52306752808 0-01 
5000 200519 A9 - -,2033142744 D-01 
8000 ‚00325 А10 = .1762226828 D-02 
10000 200515 А11 - -.4934822700 0-04 


11. Core 5 (Wavy-fin plate-fin surface 11.5-3/8W) Friction 
Factor from Kays and London (1 


Reynolds Friction 
Number Factor Polynomial Coefficients 
400 nce Al = -.1703315790 D+05 
500 2125 AZ = 71290411355 D+05 
700 1093 АЗ = -,3694551557 D+04 
900 . 088 А4 = .4359291610 D+03 
1250 . 078 А5 = -.1981584997 10-01 
1500 207/0 Аб = -.3594302044 D+01 
2000 . 0525 А? » .7340606854 D-01 
3000 ‚0525 48 = .3091640055 0-01 
4000 .0475 А9 = -.1603543052 0-02 
5000 . 043 А10 9 -.7179253410 D-04 
7000 ‚0375 All = .5035665432 р-05 


12. Core 5 (Wavy-fin plate-fin surface 11.5-3/8V) Coburn-j 
Factor from Kays and London (1) 
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Reynolds Coburn- j 
Number Factor Polynomial Coefficients 
400 .019 А] = .1845757550 р-05 
500 ‚018 A2 = -.1406316423 D+05 
700 -0175 АЗ = 4071246982 D404 
900 ‚01625 A4 2 -.4950898847 D-403 
1250 ЕЛІ» А5 = .5485777432 D+01 
1500 ‚014 Аб = ‚38283212386 10401 
2000 . 013 А7 = -.1208492523 D+00 
3000 ПОЛ A8 2 -.2934686653 D-01 
4000 0096 А9 = .1796872494 D-02 
5000 2009 А10 - .5114973685 р-04 
7000 ‚0079 А11 = -.4408556059 2-05 
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